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Inlet and Exit-Header Shapes for 
Uniform Flow Through a Resistance 
Parallel to the Main Stream 

An analytical and experimental study of flow in headers with a resistance parallel to the 
turbulent and incompressible main stream has been made. The purpose was to shape 
the inlet and exit headers , which had a large length- to -height ratio , so that the fluid would 
pass through the resistance uniformly. Analytical wall shapes and estimated total 
pressure drop through the headers were compared with experimental results. Good 
agreement between analysis and experiment was found for the cases compared. 


I here has been a great interest in the flow of fluids 
through resistances parallel to the main fluid stream, Figs. 1 and 
2. There are man3 r applications of this type of flow in such equip- 
ment as air driers and heat exchangers. The advantage of this 
flow geometry as compared with the case where the resistance is 
normal to the main stream is that a larger sized resistance can be 
used giving a lower total-pressure drop for the same mass flow 
and frontal area or, if the flow-resistance size is constant, smaller 
frontal area and less fluid holdup in the system are obtained. 
The fluid holdup refers to the amount of fluid in the line during 
steady running conditions. 

A disadvantage commonly found in this type of folded-flow 
system is that the fluid does not pass through the resistance uni- 
formly. In many applications uniform flow through the re- 
sistance is desirable primarily because then full use can be made of 
the active surfaces of the equipment, and usually the total-pres- 
sure drop through the system will be less. Several approaches to 
the problem of obtaining uniform flow exist. One method is that 
of using turning vanes [1] ; l however, these are difficult to fabricate 
and install. Because of flow separation along the vanes they are 
not always successful. Another method might be to vary the re- 
sistance along the channel; however, this may be impractical in 

1 Numbers in brackets designate References at end of paper. 
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some applications. In this study uniform flow through a porous 
resistance was obtained by shaping the inlet and exit headers. 
The fluid was assumed turbulent and incompressible and the 
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Introduction 


Nomenclature' 


a 


B 

D 

F 


f 

L 

£ 

n 


nondimensional co-ordinate nor- p 
mal to resistance for exit header, 
a /a 4 - u 

integration constant, (a t /t/*) 2 — D 
integration constant, (2 — n )/ 

(1 + 0 - n) 

, . . L (0.316) 

friction constant, 7 7 v ' 

4a t / ^2a t y- 25 

friction factor, 8 T w /pv' 2 

length of flow resistance 2/ 

defined by equation (26) 

parameter to fix exit-header shape, 

0 < n < 1 13 


nondimensional static pressure, 
p'/ipvi 2 / 2) 

nondimensional velocity in inlet 
header, u'/Vi 

nondimensional velocity in exit 
header, v'/vi 

normal velocity through resistance 

nondimensional co-ordinate along 
resistance, x' /L 

nondimensional co-ordinate nor- 
mal to resistance for inlet 
header, y' !a { 

term defined in equation (13) 


v = kinematic viscosity 
p = density 
t w = shear stress at wall 
X = defined by equation (21) 

Subscripts 

f = term to be added to account for fric- 
tion 

i = term evaluated at outlet of exit 
header or entrance of inlet header 
0 = cases where friction is neglected 

Superscript 

' = (prime) dimensional term 
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length of the header was considered large compared with the 
height. The case of shaping only the inlet header to achieve uni- 
form flow and allowing the fluid to exhaust directly on leaving the 
resistance into the atmosphere has been treated previously by 
Loeffier and Perlmutter [2]. The analysis in this paper for the 
inlet header makes similar assumptions as those used in [2]. 

The case where the fluid enters from the atmosphere through 
a resistance into an exit header which exhausts the fluid parallel 
to the resistance has been studied by Taylor [3] using both 
one and two-dimensional analyses while ignoring friction. Good 
agreement for both analyses with experiments was obtained. 
In the analysis of the exit header in this paper the one-dimensional 
analysis and similar assumptions as Taylor’s are used. It can be 
shown that it is impossible to obtain uniform flow from the at- 
mosphere through a resistance into an exit header which exhausts 
the fluid parallel to the resistance by shaping the exit header. 
This is because the exit header requires a pressure decrease in the 
direction of the outlet of the header to force out the fluid. Since 
the pressure on the upstream side of the resistance is constant, this 
causes the largest pressure difference, and hence the greatest flow 
across the resistance to be near the outlet of the exit header. 

Cichelli and Boucher [4] studied inlet and exit-header shapes 
to achieve uniform flow through the resistance. They assumed 
that the fluid leaving the inlet header carried no forward com- 
ponent of velocity. This did not prove to be true in the present 
case, as shown by experiment. They also ignored frictional ef- 
fects. Heyda [5] carried out an analytical solution for shaping 
inlet and exit headers to obtain uniform flow through the re- 
sistance using a two-dimensional approach for the exit header. 
Similar assumptions to the present paper were made except fric- 
tion was neglected and the solution carried out for linear shaped 
outer walls of the exit headers. 

In this paper the headers are designed so that the static-pres- 
sure difference across the resistance will be constant for any point 
along the resistance. Since the resistance is constant, this will in- 
sure uniform flow through the resistance at every point. This 
assumes that the streamlines are identical through the resistance. 
The problem was broken up into two separate parts, the inlet 
header and the exit header. The two problems are related by the 
fact that the pressure distribution in the exit header is the same 
as in the inlet header except for a constant difference across the 
resistance. The analysis was finally checked by experiment, Fig. 
1, and the theoretical wall shapes were compared with experi- 
mental wall shapes. Also, estimates for all-over total-pres- 
sure drops were compared with experimental over-all total- 
pressure drops. 

Analysis 

The analytical model is shown in Fig. 2. For the exit header 
for fairly large resistance, the fluid enters with no x-component of 
velocity. For very thin resistances with large open areas this 
will not be true. The one-dimensional momentum balance for 
the exit header is 


-a'dp' - 2 tJx' = pd(aV 2 ) (1) 


expected to hold exactly for this case. However, the frictional 
effects are only minor compared to the change in momentum and 
the inclusion of the friction is useful to indicate the trend and the 
magnitude of its effects. It will be shown later that the fluid 
in the inlet header enters the resistance with almost its full 
forward velocity. This would cause the shear at that resistance 
to be very small, hence the effect of friction should be smaller 
than indicated. Nondimensionalizing as follows: 


v' 

v = — , a 


o' _ p' 
a/ V (pVi 2 / 2 ) * 


and x — 


x ' 
L 


Letting F = (/ l /4)(L/ai) where 


/<- 


0.316 

(*r 


equation (2) becomes 



Fv 2 
(va ) 1 / 4 


dx = d(av 2 ) 


(3) 


Since uniform flow through the resistance is desired, from con- 
tinuity it is required that 

va = x (4) 


Substituting equation (4) into (3) yields 



At this point the pressure can be given as a function of x and 
equation (5) can be solved for the exit-header wall shape. Also, 
the exit-header wall shape can be given and then it is possible to 
solve for the pressure. In choosing the first case care must be 
taken to insure that the pressure decreases in the direction in 
which the fluid is to flow. This is because the fluid entering the 
exit header has no velocity component in the direction of the 
outlet. In order for the fluid to flow in the direction of the outlet 
the only force available would be the pressure gradient. The pres- 
sure distribution can be chosen as follows: 


V - Vi = T, ; [1 “ z 2(1 ~ n) ] (6) 

(1 - n) 

where n is arbitrary. To keep the pressure finite at x = 0, n 
must be less than 1. Differentiating equation (6) gives 

~ = 2(n - 2)*i-** (7) 


Substituting equation (7) into equation (5) yields 

da, 2a _ a\2 — n ) _F_ 

dx x x l+2n x^ K 

Exit Header With No Friction. If friction F were neglected in equa- 
tion (8), we would obtain the solution for a 


where terms in equation (1) are defined in Fig. 2 or in the nomen- 
clature. Stress t w can be defined from the Blasius law using a hy- 
draulic diameter based on infinite parallel plates in the definition 
of the Reynolds number 

8 r w 0.316 

“ t/2a'y/‘ (2) 

Actually, since there is an addition or removal of fluids and the 
wall cross section is changing, the Blasius relationship cannot be 


a 0 = x n (9) 

where the subscript 0 refers to the cases where friction is neglected. 

To prevent a 0 from becoming infinite at x = 0, n must be 
greater than or equal to zero. Therefore n is restricted between 
the limits 0 < n < 1, where the upper limit is from the previous 
paragraph. By letting n be a given positive value between these 
limits, a variety of exit-header shapes can be obtained, Fig. 3. 

Exit Header With Friction. It seems reasonable to assume that the 
effect of friction will only change the wall shape slightly; a can 
then be written as 
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Fig. 3 Exit-header shapes for various values of it, for cases of friction 
and no friction 


(l ~ CLq H - (If 


( 10 ) 


where a 0 a f . We can consider a f to be the change in wall 
shape due to friction. Substitution of equation (10) into equa- 
tion (8), and discarding the squares of a f as negligible, results in 


"7“ + — (4 — 3ra) 
dx x 



(ID 


This can be solved to give 


Fx'/* 



( 12 ) 



Fig. 4 Wall shapes fdr inlet header in counterflow, neglecting friction 


sure difference between the two headers at each x must be a con- 
stant, Using equations (7) and (14), equation (13) gives 


?-(!+«* 

ax x 


F 

x '/\ 


(2 - n) 


(15) 


The analytical solution for this equation is not obvious. It can 
he solved numerically for the inlet-header shape. If friction is 
neglected (F = 0), equation (15) can be solved to yield 


where 


JL _ D B 

y Q 2 X 2n 


D 


_ 2 - n 
1 -f 0 — n 


and 



(16) 


Adding equation (12) to equation (9), the wall shape of the exit 
header can be found for the case including friction. These results 
are plotted in Fig. 3 as dashed lines for a given value of F. For 
different values of F the resulting a f can be found by the ratio 
a f ,2 = (F 2 /Fi)a fi i. The friction causes the wall to be farther from 
the resistance than in the case for no friction, because, to maintain 
the given pressure distribution, the drop in pressure due to the 
effect of friction must be compensated for by slowing down the 
fluid by widening the channel. This causes an increase in static 
pressure to compensate for the frictional pressure loss. 

Inlet Header. The fluid entering the inlet header can be in the 
same direction (parallel flow) or in the opposite direction (counter- 
flow) to the fluid in the exit header, Fig. 2. A momentum balance 
for either case in nondimensional terms is 


The value of ajyi has a fixed lower limit for the condition of uni- 
form flow through the resistance, otherwise y 0 2 becomes negative. 
As long as 

(»:)' 2 (rrinb;) (lf “ 

y\)~ will always be positive. y 0 is plotted against x in Fig. 4. It 
can be seen that as a *-/?/; gets large the problem reduces to that of 
the fluid exhausting into a large reservoir as given in Ref. [2]. 
This condition is approached for the case of ai/yi — 4. For this 
case the inlet-header shape does not seem to depend on the exit- 
header shape n. 

Counterflow Inlet Header With Friction. Since the effect of friction 
would change the wall shape only slightly, y can be written as 


- ^ dp - F dx - d(yu 2 ) + (I — 0)ud(yu) =0 (13) 

2 \uy\ '* 

where y — y’ /a i? u — u'/vi and the other terms are made non- 
dimensional as before. In the case of parallel flow, u is positive, 
while for the case of counterflow u is negative. The term (1 — 0) 
in equation (13) can be considered the fraction of the o*-com- 
ponent of velocity of the main-stream fluid contained in the 
fluid particles passing into the resistance. The value of 0 will be 
considered a constant of the resistance and will be discussed more 
fully later on. 

Counterflow Inlet Header With No Friction. For the case w here the 
inlet-header velocity is in the opposite direction to the exit-header 
velocity (counterflow 7 ), assuming uniform flow through the re- 
sistance, from continuity it is found that 

yu — —x (14) 

The pressure gradient found from equation (7) for the exit header 
must be the same as for the inlet header because the static-pres- 


V = 2/o + V f (17) 

where y Q ^ >y s . Substituting into equation (15) and neglecting 
squares of y f yield 


dy f 

dx 


(2 - n)3y 0 *y f 


£ -f- (1 -F £) — 


Vf 


F_ 

rV< 


(18) 


Using yo from equation (16) in equation (18), y f can be integrated 
as follows: 


y f (Z)x 2(1 +^“ n > + Bf^ 
F ' x'+e 


1 

3n - 20 - - 
4 


x 


l _* 

4 (Dx 2 <d+0-") + B ) 3 /2 -f 3J5 3/i (l + (3 - n) 

Jx~ (* + f z 2(1+s ~ n) )' /! dx] + c 


(19) 
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Fig. 5 Correction to inlet-header shapes owing to friction for the counter flow case 
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( 21 ) 


Fig. 6 Inlet-header wall shape for parallel flow neglecting friction for 
values of y t /a,- = 0.5, 1, and 2 

After neglecting friction ( F = 0), this can be integrated to 
--(1 - z) 2(1 +W = 2(2 - n) I x 1 ~ 1 "(l - xf^dx + c (24' 

yo 2 J 


In Fig. 5 the values of y f /F are plotted. These values must be 
added to y Q to obtain the wall shape when the effects of friction 
are included. 

Parallel-Flow Inlet Header With No Friction. For the case where the 
inlet-header velocity is in the same direction as the exit-header 
velocity, the momentum balance is given by equation (13). From 
continuity, since flow through the resistance is constant, 


Expanding the integrand in a series yields 
(1 _ *)2(l+0 


Vo* 


£ + 




(25) 


where 


4 - 2(2 T— 

K } [_2 -2 n 3-2 in 2!(4 - 


- \)x 2 
2 n) 


u — 


1 — x 


( 22 ) 


Substituting equations (7) and (22) into equation (13) gives 

F 


20(2/3 - 1)(20 - 2) . , 

— x* + . . . 

3!(5 - 2 n) 

( — l) K (2fi)\x K 


(20 - K)\K}{2 -2 n + K) 


.] 


(26) 


ldy + (±±fi) 1 (2-n) ^ 

y 3 dx (1 — x) y 2 (1 — x) 2 
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These results are shown in Fig. 6. 

Parallel-Flow Inlef Header With Friction. Breaking y up as in equa- 
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Fig. 7 Correction to inlet-header shapes owing to friction for parallel-flow case 


tion (17), substituting into equation (23), and ignoring terms with 
y s * yield 



+ 3 y 0 2 


(2 - n) 
(1 - x) 2 



F 

(1 — x) 1 /* 


(27) 


This can be integrated as follows: 

y f exp [3y 0 2 (2 — n)(x l ~ in )/(l — x) 2 j 
F{ 1 - x) 1+ * 


-£ 


exp[3y 0 H(2 - n)(si-*")/(l - »)«] ^ (2g) 


(1 - *) 4 




The solution can be found by numerical integration of the inte- 
gral. The results from equation (28) are plotted in Fig. 7. 


Experiment 

An experimental study was undertaken to check the results of 
the analysis and study over-all pressure drops. A view of the 
test section is shown in Fig. 1. The flow enters through the upper 
inlet header. The upper wall of the inlet header, which was 
flexible, was made of Vg-in- thick Plexiglas. It could be shaped by 
means of bolts which were spaced along the sides. The entering 
height of 2 in. and the width of 16 in. were designed to approxi- 
mate a two-dimensional model. The flow left the entrance header 
through the lower surface, which was a porous plate 32 in. long. 
One or more layers of cotton cloth could be placed on top of the 
porous plate to increase the flow resistance. The exit header con- 
sisted of a wall at a fixed distance away from the resistance. 
This corresponded to an n-value of zero in the analysis. This 
distance could be made either 1 or 4 in., corresponding to an ajy^ 
respectively, of 0.5 and 2. The fluid exhausted into the room at 
the same end that the fluid had entered (counterflow). Instru- 
mentation on the test section consisted of pressure taps located 


along the center line of the upper and lower walls. The static 
pressure measured at the wall w as shown to be the same value as 
at the resistance below by static-pressure profile measurements in 
[2]. Measurements of total-pressure drop through the porous 
plates and the cloths when the fluid flow is normal to the re- 
sistance are given in [2] also. This pressure drop was found to 
follow the relationship 

5.2A H = ^ v n '«• (29) 

2 g 

where AH is pressure drop across the resistance in inches of 
water and v n ' is the velocity normal through the resistance in feet 
per second. The values of K and m were found to be 17 and 1.9 
for a porous plate only; 470 and 1.8 for a porous plate and one 
cloth; and 435 and 1.94 for a porous plate and two cloths. Dry 
air at 125 psig was available as the operating fluid. This air was 
filtered and passed through a standard ASME orifice run to 
measure the mass-flow rate. In determining the upper wall shape 
which would yield a uniform flow through the resistance, the mass 
flow was first set at some desired rate. Next the bolts, Fig. 1, 
along the sides of the flexible wall were adjusted until a constant 
pressure difference across the resistance was obtained at all points 
down the channel. If the static-pressure difference at a point 
was too high, the fluid in the upper header was accelerated by 
lowering the wall. This caused the static pressure to drop until 
the desired pressure difference was obtained. Measurements of 
the wall shape were then made by means of a depth gage. 

The experimentally obtained uniform wall shapes are com- 
pared with the theoretical wall shapes in Fig. 8. The points seem 
to be close to ft = 0, except near the entrance where the value of 
ft is higher. The porous plate without the cloth seems to stay 
close to ft — 0, even in the upstream region. The reason for this 
behavior of ft will be discussed later. 

Fig. 9 shows the static-pressure difference across the resistance 
along the channel. For the shaped wall the static-pressure dif- 
ference has been made constant along the channel by shaping the 
wall. The pressure difference was found to remain constant along 
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the channel over the flow range studied. This flow range can be 
found from the F-range in Fig. 8. Numerical calculation will 
show that changes over a large range of velocity will have a very 
small effect on the wall shape. As shown in Fig. 9, shaping the 
wall did not give constant pressure difference across the resistance 
down the channel when a x fy x was changed from 2 to 0.5. It was 
impossible to adjust the w'all to achieve uniform flow for this 
case. This agrees with the analysis for the counterflow case, 
which claims that there is a minimum value of ai///i that is above 
0.5 but below 2, equation (16a). 

To study the variation from uniform How for some different 
simple wall shape, the upper wall was made to increase linearly, 

(y' = Vix'/L). 

Since a linear wall shape gives smaller channel heights than 
found from the analysis as the flow travels downstream, the static 
pressure drops faster in this case than for the shaped wall. This is 
the reason that the pressure difference across the resistance shown 



Fig. 8 Comparison of experimental wall shape with theoretical shapes 
For counterflow case, a,;/yry = 2 


in Fig. 9 becomes smaller down the channel. Thus most of the 
normal flow passes through the resistance near the entrance of 
the inlet header for this case. For an inlet header wall a constant 
distance away from the resistance ( y ' — yi), the opposite effect 
occurs. The pressure rises downstream so that the pressure dif- 
ference across the channel increases as the fluid travels down- 
stream. For this case most of the fluid passes through the re- 
sistance in the downstream portion of the resistance. For the 
parallel flow case a similar deduction can be made. If, instead of 
the shaped wall, the height is made to rise linearly {y f — y t [ 1 — 
(x'/L)] J, then the wall will be higher than the theoretical shape 
for uniform flow, and the pressure in the inlet header will rise as 
the fluid travels downstream. This causes the largest pressure 
drop across the resistance and therefore also the largest normal 
velocity to occur close to the outlet of the exit header. 

The ratio of normal velocity through the screen at x' = 0 over 
the normal velocity at x' = L, ( v n , x ^ /v n> is approximately 

equal to the square root of the ratio of the pressure difference 
across the resistance at x f — 0 over the pressure difference at x' 
— L. This ratio was found to be only slightly dependent on flow 
rate for the straight- wall case. However, the ratio is strongly de- 
pendent on the flow resistance. The higher the resistance, the 
closer the ratio of lv n , x =\ approaches 1 . Experimentally 

for counterflow with the upper wall straight y' = y^x'/L), it 
was found that for the porous plate v n , x -o' /v Ui x _i' = 0.6. For 
the porous plate plus a cloth v nt X =</Av *-/ = 0.75. For the 
porous plate plus two cloths v n , x -VAv *-/ = 0.8. 

Total-Pressure Drop Through Test Section. A complete discussion on 
the causes of over-all pressure drops in this type of flow is given 
in [6] . Plotted in Figs. 10 , 1 1 , and 12 are the total-pressure drops 
versus flow rates for the different resistances. Shown as solid 
lines are the theoretical total-pressure losses of Kuchemann and 
Weber [7]. Their maximum estimate I assumes that the total 
head loss represents the head loss through the resistance based 
on the normal velocity v n ' as given by equation (29) plus the 
velocity head drop (p/2)(w/ 2 tf n /2 ), in turning into the screen 


RESISTANCE 



V 

Fig. 9 Pressure drop across resistance along channel 
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Fig. 10 Total-pressure drop through system for porous plate 



Fig. 1 1 Total-pressure drop through system for porous plate plus one 
cloth 



Fig. 12 Total-pressure drop through system for porous plate plus two 
cloths 

which is not considered recovered. The reacceleration in the exit 
header for this case is assumed to be achieved without losses. 
Their minimum estimate II represents either the first or second 
term of estimate I, whichever is larger. Notice that the total- 
pressure drop for the straight upper wall is higher than for the 
shaped wall. This can be explained as follows: 

For counterflow the total-pressure drop through the system can 
be considered equal to the pressure drop across the resistance at 


x r /L = 1 plus the difference in velocity head between ttie incom- 
ing and outgoing velocity at this same point. Since for the 
straight-wall case, y' — y^x'/L), the pressure difference across 
the resistance at x* /L = 1 is a maximum (as discussed earlier), 
which is higher than for the shaped w r all for the same mass-flow 
rate, a larger over-all pressure drop would be expected. A similar 
result was found in [2]. It was concluded previously that, if the 
upper wall were straight so that y' — the minimum static- 
pressure drop across the resistance would be at x'/L = 1. This 
would imply that for the same mass-flow rate the total-pressure 
drop through the system will be smaller for this case than for 
the shaped wall. Similar reasoning can be used in the parallel- 
flow case. 

Also shown in Fig. 10 is the total-pressure drop for the shaped 
wall where ajyi — 0.5. Here the total-pressure drop is much 
larger than for the other cases, since all the flow passed through 
only a small section of the resistance and thus gave a very high 
pressure drop through the resistance. 

Notice also that the case for uniform flow without an exit 
header from [2] has a higher total-pressure drop than for the case 
with the exit header. This is because the velocity head out of the 
exit header is not recovered for the case with no exit header. 

Discussion of /3. The problem of predicting the value of /3 is dif- 
ficult. One of the most complete discussions of this parameter 
appears in a paper by Soucek and Zelnick [8] , and in the discus- 
sion at the end of that paper by McNown. They came to the 
conclusion that (3 was dependent on the surface characteristics of 
the resistance and on the ratio of the exiting velocity to the main- 
stream velocity. McNown felt that /3 could never be greater 
than 0.5. Cichelli and Boucher [4] in their analysis assumed a 
value of (3 of 1. Heyda [5] assumes a jS of 0. 1 — /3 in equation 
(13) can be considered as the fraction of the forward velocity of 
the main-stream fluid contained in the fluid leaving the inlet 
header. 0 would be 0 if the exiting fluid had an ar-component of 
velocity equal to that of the main stream. {3 would be 1 if the 
fluid particle turned and left the inlet header with no forward 
component of velocity. On the basis of smoke streamline studies 
made by Comtois [1], it was found that the fluid enters the screen 
with very little turning. This would be expected, since no pres- 
sure gradient was found normal to the resistance. Thus there 
would be no force to turn the particle until it enters the re- 
sistance. This means that the particle leaves the inlet header 
with its full forward momentum. The experimental results ob- 
tained here show a (3 greater than zero near the entrance of the 
inlet header which drops to 0 farther downstream. The reason for 
this is that initially the boundary layer is being removed. The 
average velocity across the boundary layer is less than the average 
velocity across the header. Therefore, although the boundary- 
layer fluid leaves with its full forward component of velocity this 
is less than that of the mean velocity across the channel. How- 
ever, as this boundary layer is removed, approaches zero. 

Conclusions 

A study was made of turbulent incompressible flow through re- 
sistance parallel to the stream direction for the purpose of obtain- 
ing uniform flow through the resistance. By shaping the inlet 
and exit headers, which had a large length-to-height ratio, that 
was found to be possible. Theoretical and experimental wall 
shapes for uniform flow are given and show good agreement. 
These wall shapes were found to stay constant over a large veloc- 
ity range. Experimental results of the effect of using some simple 
wall shapes different from those calculated theoretically are also 
shown. 

Measurements of over-all total-pressure drops were taken and 
compared with theoretical estimates with good agreement. It 
was shown that for certain shapes giving nonuniform flow it was 
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possible to obtain lower over-all pressure drops than for the case 
of uniform flow through the resistance. 
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DISCUSSION 
J. F. Heyda 2 and C. D. Fulton 2 

We are in substantial agreement with the author but offer the 
following remarks. 

There is no reason to believe in the existence of the quantity 
beta, which implies in Equation (13) that if any fluid becomes re- 
tarded its energy will be added to the remaining fluid. This is 
impossible in either potential or viscous flow. 

Beta may have originated in the notion that the fluid would 
start turning while still in the inlet header because the turn must 
be complete immediately afterward and in subsonic flow nothing 
should happen abruptly. In potential flow, even if there were 
such a turning, there could be no transmission of energy from one 
stream tube to another, and therefore beta could not exist. But 
more than that, there is no such turning. A map of curvilinear 
squares — a Laplacian field — constructed in the inlet header proves 
this. The porous resistance dictates only the position at which 
each streamline will be swallowed and not the angle at which it 
will approach. The angle is, in fact, equal to that of a line drawn 
from the far end of the resistance to the header wall opposite the 
point of swallowing. Thus the flow strikes the resistance with full 
force and delivers its momentum to it after leaving the inlet 
header. 

If the surface of the resistance contains large holes, vanes, or 
other configurations, the potential flow pattern will show repeti- 
tive stagnations and regions of increased velocity. But the 
average velocity just outside the holes or vanes wall be the same 
as in the inlet header, and so will the average pressure. Thus, 
whether the microscopic or the macroscopic view be taken, the 
results will be the same in the inlet header in potential flow 7 , and 
beta will not exist. 

Nor can beta exist because of frictional effects. If viscous 
shear retards the adjacent fluid in the header before it is swal- 

2 Aircraft Nuclear Propulsion Department, General Electric Com- 
pany, Cincinnati, Ohio. 


lowed, the momentum extracted from this fluid will again appear 
as a forward force in the solid material of the resistance. The 
total force on the material will be the same with or without 
viscous shear, since the material is bound to receive all the momen- 
tum one way or another. The momentum extracted by shear 
will not appear as a forward push on the main stream of the 
header. 

In the turbulent flow, viscosity can also create eddies of re- 
tarded fluid which will attempt to propagate into the main stream 
at an acute angle. In some cases, that angle will be smaller than 
the angle at which the flow is leaning into the resistance, and all 
the eddies will be swallowed. In other cases, a few eddies will 
propagate and will thereby reduce the pressure in the header. 
They can be accounted for by adding a small increment to the 
value of F. What Equation (13) actually states, instead, is that 
if friction underlies the existence of beta, then friction will pro- 
duce a rise in pressure, which is absurd. 

The notation used by the author in defining dimensional and 
nondimensional variables is not quite consistent in the Nomen- 
clature. This produces difficulties later. Thus, to be consistent, 
di should be replaced by a/, and aq would not need to be used since 
its value w 7 ould then be unity. Similar remarks apply to and y/. 

In the exit header, two types of flow may be conceived. One is 
mixed or one-dimensional and can be generated if the holes in the 
plate squirt out very energetic jets or if some other strong stirring 
mechanism is at work. This corresponds to Equation (3). 

The other type of flow is stratified and two-dimensional. It 
can occur if the flow 7 leaves the resistance in a relatively quiescent 
state. Each stream tube retains its identity and moves in a 
curved, accelerating path toward the outlet. The velocity profile 
is steeply tilted. We have analyzed this flow 7 using a one-dimen- 
sional pressure field, w 7 hich is justifiable both analytically and 
experimentally. The result is the following equation, here written 
in dimensional form but omitting the author’s prime notation. 

a(x) = v n Vp/2 j , dt ~ =• 

Joy/ p(t) — p(x) 

The integral computes the header width by finding the vertical 
widths of stream tubes originating at positions t when they reach 
position x. Standing at one position x, the variable t is run from 
0 to x and the integral computed. Then x is changed and the 
process repeated. This is a Volterra improper integral. As t ap- 
proaches x, the vanishing of the denominator reflects the fact 
that each stream tube starts out as a second-degree parabola that 
takes up a great deal of space. This voluminousness is partly 
compensated for by the narrowness of the high-speed tubes near 
the header wall. The equation can be converted into Abel’s in- 
tegral equation. It can also be inverted approximately for p 
by tne methods of Reference [5] . It can be elaborated for com- 
pressibility and variable flow rate through the resistance. It can 
be solved in closed form for certain simple functions. 

The results of the equations for mixed and stratified flows in 
the exit header are not greatly different although the equations 
appear very different and the mathematics are more difficult in the 
improper integral. For a given header shape, the form of the 
pressure curve is the same and the pressure gradient is merely 
somewhat steeper in the stratified flow. For this reason, one 
may elect to use the mixed flow equation corrected by a multiplier 
to handle stratified flow, or any flow between the two. Many ex- 
perimental results agree closely with stratified flow and show the 
associated velocity profile. The following table shows how the 
two flows compare. 
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Exit header shape 
a ~ x Triangle 
a ~ \/ x Parabola 
a = C Rectangle 


Shape of pressure curve 
p( 0) - p(x) 

C — log 1/x Logarithmic* 


x Linear 

x 2 Square law 


End-to-end pressure fall in 
velocity heads at mean outlet 
velocity 

Mixed Stratified 

flow flow 

4.60* 7.23* 

3 4 

2 2.47 


* In the triangular header the pressure rises to an infinite peak at x = 0. The constant C 
is indefinitely large. The tabulated pressure falls are for the latter 90 per cent of the header 
length, for this case only. 


The shapes of the inlet and exit headers can be related directly 
by eliminating the pressure. For a given exit header, there is a 
multiplicity of inlet headers nested one inside another. The larg- 
est possible inlet header in counterflow is the one producing 
stagnation at its closed end. It requires the least possible inlet 
energy and may be called the economic shape. Omitting friction, 
that shape is given by the following equations when mixed flow 
occurs in the exit header. One equation is the inversion of the 
other. The choice of which one to use depends upon which header 
is given at the beginning. These equations may be regarded as 
either dimensional or nondimensional. 

fe]’ = Li]' + 2 £[i>]'f 


When stratified flow occurs in the exit header, the correspond- 
ing equation is the following, which is not easy to invert. There- 
fore one must start with y{x). There is a multiplicity of shapes 
y(x) that give the same a(x). Economic shapes are those having 
appreciable width at the closed end — i.e., an infinite slope there. 


a(x) = 



dt 


V x i 2 _ r _tjy 
.y(x) J 


The discusser has performed experiments 4 on flow configura 
tions similar to those investigated by the author. However, the 
transverse resistances were long slots which permitted the fluid 
entering the exit header to do so with the full (inlet header) x- 
component of velocity. In fact, as far as the inlet header is con- 
cerned, it was found that the factor /3 was slightly negative. 
These conditions represent the opposite case to the author’s and 
it is interesting to see what conclusions can be drawn using the 
author’s method of analysis. 

Assuming that the fluid enters the exit header with an x-com- 
ponent of velocity u', the nondimensional momentum equation 
(for parallel flow) corresponding to the author’s equation (3) 
becomes with his symbols, 

d Fv 2 

dp . dx + u d(av) = d(av 2 ) (3)' 

2 (va) /4 

The continuity equations for parallel flow and uniform transverse 
discharge become 

va — x (4)' 

uy = 1 — x (22)' 

Substituting the continuity equations into equation (3)', 



The results of the foregoing three equations are shown in the 
following table: 


Exit header 
shape 
Triangle 
Parabola 
Rectangle 


Economic inlet 
header shape 
Impossible 
Parabola 
Rectangle 


Economic inlet header width 
Exit header width 
Mixed Stratified 

exit flow exit flow 

(L577 (L5 

0.707 0.636 


J. F. Thorpe 3 


The author has investigated a very interesting problem and 
presents information which is valuable to the fluid mechanics 
engineer. One particularly useful aspect of this investigation is 
that it defines a class of flow problems to which the one-dimen- 
sional approximation can be applied with acceptable results. 

The author has assumed that the fluid enters the exit header 
with no x-component of velocity. This assumption is valid for 
the high resistance of the porous cloths which he used. As a 
result of this assumption he finds that it is impossible to main- 
tain a constant pressure in the exit header. It was also found 
that the momentum factor /?, which is applied to the fluid leaving 
the inlet header, was essentially zero or slightly positive for the 
resistances tested. 


3 Senior Engineer, Reactor Development and Analysis Department, 

Bettis Atomic Power Laboratory, Westinghouse Electric Corpora- 
tion, Pittsburgh, Pa. Assoc. Mem. ASME. 


Now if friction is neglected ( F = 0), it can be seen that equation 
(5)' admits of a solution in the case n = 1. That is if 


a = a 0 = x n = x; ( n ~ 1), 


Then 


— dp + — — dx - dx 

2 y 


(9)' 


All that is necessary is to take 


Then 


y = y 0 = 1 - x 


dp = 0; u — v — 1 


(30) 

(31) 


Thus it is possible to keep the static pressure constant in the exit 
header by shaping both the inlet and exit headers with straight 
walls. 

If the solution outlined above is to be entirely consistent then 
equation (13) must reduce to dp = 0. Otherwise the assumed 
uniform transverse flow will not occur. With F = = 0 it is 

seen that equation (13) is satisfied with dp — 0 (recall u — 1). 

The discusser has measured the pressure distributions 4 for the 
flow configuration defined by equations (9)' and (30). It was 
found that the pressure distributions in the inlet and exit headers 


4 J. F. Thorpe, “A Parallel Duct Flow Problem,” Doctoral Dis- 
sertation, University of Pittsburgh, 1960. 
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were constants. From this result, it was concluded (as did the 
present author) that the effect of friction has little to do with the 
wall shape for uniform transverse discharge. 

In the experiments performed by the discusser, the transverse 
openings were slots in a vicinity where the local fluid velocity 
was greater than the average velocity in the inlet header. From 
the measured velocity profiles the factor 0 was estimated to be 
about —0.035. In calculating the axial pressures it was found 
that best results were obtained if this negative value was actually 
used in the calculations. Thus, it would appear that the pres- 
ent method of handling this momentum factor is valid. 


It is desired to minimize ajyi so as to reduce p to tai. 
Rewriting ( 16) as follows for 0 = 0 


Vo 2 


M(r^)(*' ,, -"- , ) + (S)'] (3 ‘> 


Since 1 /t/o 2 can never be negative the term in the bracket [ ] 

must be positive for all x. Since it is most likely to be negative 
at x = 0 the result for the bracket yields 


(*)' > (r^D 


(35) 


Author’s Closure 

The discussers' comments were very worth while and have 
added a great deal to our understanding of the present problem. 
Some points have been raised which should be enlarged upon. 
A 0 not equal to zero does not always imply that some of the 
energy of the leaving fluid is added to the remaining fluid in the 
inlet header. If there is some nonuniform entering velocity 
profile in the inlet header, the mean velocity near the porous 
wall will not be equal to the mean velocity across the channel 
and when this fluid near the wall is removed with all its for- 
ward momentum, 0 will not be zero. 

The discussers,’ conclusions on the shapes with the lowest total 
pressure drops were very interesting. For the author’s case this 
economical shape can be arrived at as follows. The total pres- 
sure drop through the system would be in nondimensional form 


Ptotai piulet Poutlet “b 



(32) 


The smallest possible value of a,-/^ then will occur at an n of 0. 
Then ^/a t - = 0.707 and both the inlet header and exit header 
will be rectangles as found by the discussers. The total pressure 
drop w ill then be 

Ptotai = R + 1 (36) 

Since from (33) the smallest value of possible is desired 
for the most economical shape this would occur in mixed flow 
with rectangular headers, as shown in the discussion. 

For the case of parallel flow from (32) and (6) 

2 — 71 ( CL; \ 2 

P total = R + : + ( — ) - 1 (37) 

1 ~ n \iji / 

This will be a minimum for the rectangular exit header (n = 0) 
then 

P total = R -f- 1 -{- (38) 


For the counterflow case p inlet — pout let will differ by some 
constant of the screen resistances R. By using continuity (32) 
results in 



The most economical shape would be the largest inlet and 
smallest outlet feasible. However, only in the limit (a^/i/i — ► 0) 
will it be as economical as the counterflow 7 case (36). 

Finally there is some question whether 0 actually equals 
—0.035 for the slotted resistance case. This small effect could be 
possibly due to friction in the header. 


Printed in U. S. A. 


370 / SEPTEMBER 1961 


Transactions of the AS ME 


